Charm elliptic flow at RHIC 
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Charm elliptic flow in heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) is studied 
in a multiphase transport model. Assuming that the cross section for charm quark scattering with 
other light quarks is the same as that between light quarks, we find that both charm and light quark 
elliptic flows are sensitive to the value of the cross section. Compared to that of light quarks, the 
elliptic flow of charm quarks is smaller at low transverse momentum but approaches comparable 
values at high transverse momentum. Similar features are seen in the elliptic flow of charmed 
mesons as well as that of the electrons from their semileptonic decays when the charmed mesons 
are produced from quark coalescence during hadronization of the partonic matter. To describe the 
large electron elliptic flow observed in available experimental data requires a charm quark scattering 
cross section that is much larger than that given by the perturbative QCD. 

PACS numbers: 25.75.Ld, 24.10.Lx 



I. INTRODUCTION 

Because of their large masses, heavy quarks are pro- 
duced very early in ultra-relativistic heavy ion collisions 
through hard collisions between nucleons. Their initial 
momentum spectra can thus be described by the pertur- 
bative quantum chromodynamics (pQCD). How their fi- 
nal spectra deviate from the initial ones depends on their 
interactions in the initial partonic matter and the mech- 
anism that converts them to hadrons as well as their sub- 
sequent interactions in the hadronic matter. Using two 
extreme scenarios for the charmed meson spectrum, i.e., 
pQCD without final-state interactions and completely 
thermalized hydrodynamics with transverse flow veloc- 
ity field, it was, however, found in Ref.pJ that the trans- 
verse momentum spectra of the electrons from their decay 
are both consistent with that measured in Au+Au colli- 
sions at RHIC J31 ■ The two scenarios give, on the other 
hand, very different elliptic flows for charmed mesons as 
well as their decay electrons, when the hadronization 
of charm quarks is modeled by coalescence with light 
quarks jjjj. Study of heavy flavor particle elliptic flow 
in ultra-relativistic heavy ion collisions thus provides the 
possibility of probing the interactions of heavy quarks in 
the partonic matter formed in these collisions and their 
hadronization mechanism. 

The study in Rcf. 3] is schematic as the charm quark 
elliptic flow in the scenario of complete thermalization is 
assumed to be the same as the fitted light quark elliptic 
flow from experiments. To understand quantitatively the 
effect of final-state interactions on the charmed meson 
spectrum and elliptic flow, we use in the present study a 
multiphase transport (AMPT) model that includes both 
initial partonic and final hadronic interactions as well as 
the transition between these two phases of matter. As- 
suming that the cross section for charm quark scatter- 
ing with light quarks is the same as that between light 



quarks, we find that the charm quark elliptic flow is sen- 
sitive to the value of the cross section and exhibits a 
transverse momentum dependence that is very different 
from that of light quarks. While the elliptic flow of charm 
quarks at low transverse momentum is much smaller than 
that of light quarks as expected from the mass ordering 
of particle elliptic flows predicted in the hydrodynamic 
model, the two become comparable at high transverse 
momentum. For final charmed hadrons and their decay 
electrons, their elliptic flows are found to follow closely 
that of charm quarks and thus carry information on the 
charm quark interactions in the initial partonic matter. 
Our results are similar to that found in Ref.jj] based on 
Molnar's Parton Cascade (MPC) model. On the other 
hand, the value of the charmed meson elliptic flow from 
our study is much larger than that found in Ref. [j| based 
on the Hadron-String Dynamics (HSD) model, in which 
the initial dense matter consists of strings and hadrons in- 
stead of partons as in the AMPT and MPC models. The 
thermalization and elliptic flows of charm quarks in a hy- 
drodynamically expanding quark-gluon plasma have also 
been studied in a Langevin model, and they are found to 
be sensitive to the charm quark diffusion coefficient [||. 
We note that preliminary results from the present study 
have been reported in Ref.Q. 



This paper is organized as follows. In Sect. [HJ we 
give a brief review of the AMPT model and then extend 
it to include the dynamics of charm particles. Results 
from the AMPT model on both the transverse momen- 
tum spectra and elliptic flows of charm particles are pre- 
sented in Sect. IIIII Finally, a summary is given in Sect. 
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II. THE AMPT MODEL 
A. A brief review 

The AMPT model H, 0, E3, El E3, O is a hybrid 
model that uses minijet partons from hard processes and 
strings from soft processes in the heavy ion jet interaction 
generator (HIJING) model as the initial conditions 
for modeling heavy-ion collisions at ultra-relativistic en- 
ergies. Since the initial energy density in Au+Au colli- 
sions at RHIC is much larger than the critical energy 
density at which the hadronic matter to quark-gluon 
plasma transition would occur 0, fl3l fl5| , we use the 
version which allows the melting of initial excited strings 
into partons |l6[ . In this version, hadrons that would 
have been produced from the HIJING model are con- 
verted to valence quarks and/or antiquarks. Interactions 
among these partons are described by Zhang's parton 
cascade (ZPC) model ^J- At present, this model in- 
cludes only parton-parton elastic scatterings with an in- 
medium cross section given by the perturbative QCD, 
i.e., 

da p 9na 2 s f o?\ 1 

where a s = 0.47 is the strong coupling constant, and s 
and t arc the usual Mandelstam variables. The effective 
screening mass /x depends on the temperature and den- 
sity of the partonic matter but is taken as a parameter 
in ZPC for fixing the magnitude and angular distribu- 
tion of parton scattering cross section. Since no inelastic 
scatterings are included in the ZPC model, only quarks 
and antiquarks from the melted strings are present in the 
partonic matter. The species independence of the above 
cross section compensates, however, for the absence of 
gluons in the early stage. 

After partons stop scattering, they are converted to 
hadrons using a coordinate space quark coalescence 
model, i.e., two nearest quark and antiquark are com- 
bined into mesons and three nearest quarks or antiquarks 
are combined into baryons or anti-baryons that are clos- 
est to the invariant masses of these parton combinations. 
This coalescence model is somewhat different from the 
ones that are based on the overlap of the hadron quark 
wave functions with the quark distribution functions in 
the partonic matter and used extensively for studying 
the production of hadrons with intermediate transverse 
momenta 0, 0, |2(j ■ The final-state scatterings of pro- 
duced hadrons in the AMPT model are described by a 
relativistic transport (ART) model [2l|. 

Using parton scattering cross sections of 6-10 mb, the 
AMPT model with string melting is able to reproduce the 
centrality and transverse momentum (below 2 GeV/c) 
dependence of hadron elliptic flows ^(| and higher-order 
anisotropic flows |22| as well as the pion interferometry 
[23| measured in Au+Au collisions at yfs = 130 AGeV 

at RHIC millll^. Tt has also been used for stud y m S 



the kaon interferometry [2Jj in these collisions as well as 
many other observables at y/s = 200 AGeV pll^ . 



B. Including charm quarks and charmed mesons 

To include charm particles in the AMPT model, we 
first generate the initial charm quark distributions using 
information obtained by the STAR collaboration from 
d+Au collisions at y/J^N = 200 GeV Since 
final-state interactions are negligible in this collision, the 
measured D meson spectrum is related via the binary col- 
lision scaling to that from the nucleon-nucleon collision at 
same energy. In panel (a) of Fig. ^ we show the extracted 
transverse momentum spectrum of D° mesons in the ra- 
pidity interval \y\ < 1 that are reconstructed from either 
the Kit invariant mass (solid squares) or the Kirp invari- 
ant mass (open circles). Following Refs. [M 123. IjfflL l3ll| . 
we parameterize this spectrum by a power law, 



2Trp T dp T dy \ 3.25(GeV/c) J ' w 

as shown by the solid curve in Fig. ^a). The positions 
of these charmed mesons are then distributed according 
to those of initial binary nucleon-nucleon collisions. 

As for light hadrons produced from string fragmenta- 
tion in the string melting scenario of the AMPT model, 
we dissociate these D mesons into charm and anti-charm 
quarks according to their valence structures to obtain the 
initial charm quark distributions. The charm quarks then 
propagate along straight-line trajectories from the posi- 
tions where they are produced in nucleon-nucleon colli- 
sions for a duration given by a formation time that is 
taken to be the inverse of the D meson transverse mo- 
mentum. 

Since the cross section for charmed meson production 
in a nucleon-nucleon collision extracted from the STAR 
d+Au data is about 1.3 mb and is about a factor of 30 
smaller than the nucleon-nucleon inelastic cross section 
(~ 40 mb), the number of charm particles produced in 
relativistic heavy ion collisions is small. To increase their 
statistics in the Monte Carlo treatment of the AMPT 
model, we apply the perturbative method that have been 
extensively used in transport models for studying rare 
particle production in heavy ion collisions at lower en- 
ergies (33, |33|| . Specifically, we scale the charm produc- 
tion cross section by a factor that equals the ratio of the 
nucleon-nucleon inelastic cross section to that for charm 
production. Each charm particle thus carries a proba- 
bility that equals the inverse of this scale factor. As ex- 
perimental data indicate that charm production in heavy 
ion collisions at RHIC is consistent with the binary colli- 
sion scaling, we only allow charm quarks to be produced 
in initial nucleon-nucleon collisions. Rescattering of the 
charm quarks in the initial partonic matter and that of 
resulting charmed mesons in later hadronic matter is, 
however, included through the AMPT model. Although 
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momenta of these charm particles are modified after scat- 
tering with other more abundant particles, their effects 
on the latter are neglected due to the smaller number of 
charm particles in each event. 

Both the bulk dynamics of initial partonic matter and 
the time evolution of charm quarks are then described 
by the ZPC model using current quark masses, i.e., 9.9 
MeV for d quark, 5.6 MeV for u quark, 199 MeV for 
s quark, and 1.35 GeV for c quark. In this study, we 
use two values for the cross section: <r p = 3 mb esti- 
mated from the perturbative QCD and o v = 10 mb for a 
strongly interacting QCD. Charm rescattering cross sec- 
tions have been studied before in 0, |3f| , and the 3 mb 
cross section gives an upper bound for the correspond- 
ing energy dependent charm rescattering cross sections 
when singularities are regulated by the screening mass. 
We will use, however, the energy independent cross sec- 
tion to study whether experimental data require a cross 
section much larger than the perturbative one. The dif- 
ferential cross section then follows Eq. (JJ with s replaced 
by the maximum squared momentum transfer for charm 
rescatterings. As for the hadronization of light quarks, 
charm quarks are converted to hadrons using the coor- 
dinate space quark coalescence model described in the 
previous section. Because of the perturbative method 
used for treating charm quarks, charmed mesons acquire 
same probabilities as those of their parent charm quarks. 
Although the scattering cross sections of D mesons with 
other hadrons have been studied in an effective hadronic 
model 36] , we use in the present study instead the same 
cross section as that for parton-parton scattering. We 
have found that hadronic scattering does not affect much 
the final charmed meson spectrum and elliptic flow in 
heavy ion collisions at RHIC. To compare with mea- 
sured electrons from charmed meson decays, we include 
D mesons that are both directly produced as well as from 
the decay of D* mesons at freeze out. 



III. RESULTS 

A. Charm transverse momentum spectra 

The results from the AMPT model on the transverse 
momentum distributions of charm particles at midrapid- 
ity from minimum bias Au+Au collisions at ^/snn = 200 
GeV are shown in Fig. ^ Panel (b) shows the ini- 
tial charm quark transverse momentum spectrum (dot- 
ted curve) together with final ones obtained using charm 
quark scattering cross sections of 3 mb (dashed curve) 
and 10 mb (solid curve). The initial charm quark trans- 
verse momentum spectrum is softer than that of initial 
D° mesons shown in panel (a) , and this is due to the effec- 
tive quenching introduced by converting D mesons into 
their valence quarks. Partonic scattering leads to fur- 
ther quenching of charm quark spectrum, with the effect 
increasing with increasing charm quark scattering cross 
section. The quenching effect is, however, reversed when 
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FIG. 1: (Color online) Transverse momentum distributions at 
midrapidity in minimum bias Au+Au collisions at ^/sjvjv = 
200 GeV: (a) D° mesons per binary nucleon-nucleon collision 
extracted from d+Au collisions in the STAR experiment |3(il| 
with solid squares and open circles denoting data and the solid 
curve from the parameterization (Eq.©); (b) Initial (dotted 
curve) and final (dashed curve for a v — 3 mb and solid curve 
for <t p = 10 mb) charm quarks from AMPT; (c) Final charm 
quarks (thin dashed and solid curves) and D mesons (thick 
dashed and solid curves); and (d) Electrons from D meson de- 
cay from AMPT (dashed and solid curves) and the PHENIX 
experiment |2jj (solid circles). 



charm quarks and light quarks are recombined into D 
mesons via the coalescence model as shown in Fig.^c), 
where the charm quark spectra at the end of partonic 
phase (thin curves) and the D meson spectra at the end 
of hadronic phase (thick curves) are compared. We note 
that the final D meson number in midrapidity is about 
20% larger than that of charm quarks in the same rapid- 
ity, and they are produced from charm quarks at larger 
rapidity as a result of the coordinate space coalescence. 

Although the D meson spectrum in Au+Au collisions 
has not been measured at RHIC, that of the electrons 
from D meson semi-leptonic decay is available from the 
PHENIX collaboration |3j| as shown by solid circles in 
panel (d) of Fig. ^ They are very well described by the 
AMPT model with both charm quark scattering cross 
sections of 3 mb (dashed curve) and 10 mb (solid curve). 
Similar to that found in Ref.Q], these experimental data 
thus cannot distinguish the dynamics of charm quarks in 
the partonic matter. This is not surprising as the decay 
process significantly softens the transverse momentum 
spectrum of electrons with respect to that of D mesons. 
We note that in this comparison the charm production 
cross section in a nucleon-nucleon collision is taken to be 
600/ib in order to be consistent with the PHENIX mea- 
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surement. 



B. Charm elliptic flow 

We have also studied the elliptic flows of charm quarks 
and D mesons, i.e., the second Fourier coefficient w 2 i n 
the decomposition of their transverse momentum distri- 
butions with respect to the azimuthal angle <f> in the re- 
action plane |38j . 
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2v n (p T ,y) cos(n</>) 



(3) 

Because of the symmetry cf> ^ —<j) in the collision ge- 
ometry, no sine terms appear in the above expansion. 
The anisotropic flows v n generally depend on the parti- 
cle transverse momentum and rapidity, and for a given 
rapidity the anisotropic flows at transverse momentum 
Pt can be evaluated from 



Vn(pr) = (cos(n0)) 



(4) 



where (• • •) denotes average over the azimuthal distri- 
bution of particles with transverse momentum px- The 
elliptic flow vi can further be expressed in terms of the 
single-particle average: 



V2{pt) 




(5) 



where p x and p y are, respectively, projections of particle 
momentum in and perpendicular to the reaction plane. 
It has been shown that the elliptic flow is sensitive to the 
early dynamics of produced matter in relativistic heavy 

ion collisions [H El E3> ES S H > and lt is thus an 
especially robust observable for studying the interactions 
of early produced particles like charm quarks in the par- 
tonic matter. 

In Fig. we show the transverse momentum depen- 
dence of elliptic flows of light, strange, and charm quarks 
as well as those of D mesons and electrons from their de- 
cays in minimum bias Au+Au collisions at y^sjvjv = 200 
GeV. As shown in the upper panels (a) and (b), the i>2 
displays a clear mass dependence, with the lighter quark 
having a larger value of vi at low pt and reaching the 
maximum value at a smaller pt- This mass effect is 
particularly strong between the heavy charm quark and 
other light quarks. On the other hand, the maximum 
value of vi is similar for all quarks, exhibiting thus a 
very weak quark mass dependence. It is, however, sen- 
sitive to the parton scattering cross section, increasing 
from about 6% to about 9% when the parton scattering 
cross section changes from a v = 3 mb to 10 mb. The 
mass dependence of quark vi at lower pt is similar to 
the mass ordering of particle elliptic flows observed in 
hydrodynamic models. Our results thus imply that with 
a parton scattering cross section larger than a p = 3 both 
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FIG. 2: (Color online) Elliptic flows of light, strange, and 
charm quarks (upper panels), D mesons, and their decay 
electrons (lower panels) in minimum bias Au+Au collisions 
at ^snn = 200 GeV for parton scattering cross sections of 
Op = 3 mb (left panels) and a p = 10 mb (right panels) . 



low px charm quarks and light quarks probably approach 
thermal equilibrium in heavy ion collisions at RHIC. 

From the lower panels of Fig. |3 it is seen that the D 
meson elliptic flows (open diamonds) follow closely corre- 
sponding charm quark elliptic flows shown in the upper 
panels of Fig. although they are slightly shifted to 
higher pt and also have slightly higher values relative to 
those of charm quarks. Also shown in the lower panels of 
Fig.|3are the elliptic flows of electrons from D meson de- 
cays (solid triangles), and they also follow closely the cor- 
responding D meson ones. It is interesting to see that at 
around pt — 300 MeV, the electron elliptic flow becomes 
slightly negative, indicating that on average these elec- 
trons are moving perpendicular to their parent D mesons. 
At px above 2 GeV, the electron elliptic flow is slightly 
smaller than the D meson elliptic flow. Although the 
electron pt differential elliptic flow is almost identical to 
that of D mesons, the integrated electron elliptic flow is 
much smaller than that of D mesons as the electron trans- 
verse spectrum is much softer than that of D mesons. We 
find, e.g., v 2e = 0.55 ± 0.01% < v 2D = 2.73 ± 0.04% and 
v 2e = 0.90 ± 0.01% < v 2D = 3.91 ± 0.04% for parton 
cross sections of 3 and 10 mb, respectively. 

In Fig. |31 we compare the elliptic flow of electrons from 
D meson decay obtained from the AMPT model with 
available experimental data from the PHENIX [45j (solid 
circles) and STAR 0] (open circles) collaborations. The 
theoretical results are shown by the dashed and dotted 
curves for charm scattering cross sections of 3 and 10 
mb, respectively. Compared with the experimental data, 
the calculated charm flow has a similar magnitude at 
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FIG. 3: (Color online) Transverse momentum dependence of 
the elliptic flows of electrons from D meson decays from the 
AMPT model for different parton scattering cross sections 
of 3 (dashed curve) and 10 (dotted curve) mb with angular 
distribution given by Eq.Q as well as 10 mb with isotropic 
angular distribution (solid curve) from minimum-bias Au+Au 
collisions at ^/snn = 200 GeV. Available experimental data 
are from the PHENIX (solid circles) and STAR 46] (open 
circles) collaborations. 



1 < pr < 1.5 Gev/c but is somewhat smaller at other 
transverse momenta. 

In Ref.jj], the elliptic flow of electrons from D me- 
son decay has also been studied in the coalescence model 
based on charm and light quark elliptic flows that are 
similar to ours shown in Fig. [21 and its value at high 
transverse momentum is larger and thus closer to the 
experimental data. The smaller charm elliptic flow from 
our study is partially due to the quark spatial anisotropy, 
particularly the negative spatial eccentricity S2, that is 
present in the dynamical AMPT model but is neglected 
in Ref . ;3J . As shown in Ref . 47] , a negative quark S2 re- 
duces the elliptic flow of produced hadrons in a dynamic 
coalescence model such as the one used in the AMPT 
model. 

As mentioned in Section ^ the magnitude of parton 
scattering cross sections used in the AMPT model is de- 
termined by the screening mass or parameter /i shown 
in Eq.Q. Since the larger cross section a p — 10 mb is 
obtained by using a smaller value for /i, its angular dis- 
tribution becomes more forward peaked than in the case 
of 3 mb, thus reducing the effect of increasing cross sec- 
tion on the transport properties of quarks. On the other 
hand, it has recently been suggested that the scatter- 
ing between charm and light quarks might go through a 
quasi colorless charmed resonance, leading to a cross sec- 
tion that is not only large but also isotropic ji^. To see 
the effect of such a scattering mechanism, we have car- 



ried out a calculation with an isotropic charm-light quark 
scattering cross section of 10 mb, and the elliptic flow of 
electrons from the decay of resulting D mesons is shown 
by the solid curve in Fig. [3J It is seen that this indeed 
further enhances the charm elliptic flow but not enough 
to reproduce available experimental data. If the charm 
elliptic flow remains large when more accurate data are 
available, then the charm quark scattering cross section 
will be even larger than what have been used in present 
study. 

We note that a naive or additive quark coalescence 
model [i^. l50| . which ignores the effect of momentum 
distribution of quarks inside hadrons [fill IH^ , the non- 
negligible local spatial anisotropy > an d the large local 
directed flow |53J, has been sometimes used for calculat- 
ing the hadron elliptic flows from those of quarks. In this 
case, the elliptic flows of hadrons are then simple sums 
of the elliptic flows of their constituent quarks. Using 
this simplified quark coalescence model, the elliptic flow 
of electrons from D mesons produced from the charm 
and light quark elliptic flows shown in Fig.|21for a parton 
scattering cross section of 10 mb is comparable to that 
seen in available experimental data. 



IV. SUMMARY 

We have studied the transverse momentum spectra and 
elliptic flows of charm particles and the electrons from 
charmed meson decay in heavy ion collisions at RHIC in 
the AMPT model by treating the screening mass in the 
in-medium pQCD scattering cross section between light 
quarks as a parameter and assuming that charm quarks 
have similar cross sections. Using an initial charm quark 
transverse momentum spectrum obtained from the dis- 
sociation of the empirically determined charmed meson 
spectrum in nuclcon-nucleon collisions at same energy, 
we find that the final transverse momentum spectrum 
is softer than that in nucleon-nucleon collisions as a re- 
sult of the quenching effect due to charm quark rescat- 
tering in the partonic matter. Although the quenching 
effect increases with increasing charm quark scattering 
cross section, the transverse momentum spectrum of elec- 
trons from resulting D meson decays does not depend 
much on the parton scattering cross sections as a result 
of the stronger quenching effect due to the decay pro- 
cess, and results using both cross sections of 3 and 10 mb 
are consistent with available experimental data from the 
PHENIX experiment. 

The elliptic flow of charm quarks is, however, sensi- 
tive to the parton cross section as that of light quarks. 
With respect to the light quark elliptic flow, the charm 
quark elliptic flow is smaller at low transverse momen- 
tum but reaches a similar value at high transverse mo- 
mentum. With the coordinate space quark coalescence 
model for hadronization in the AMPT model, the elliptic 
flow of D mesons is found to follow that of charm quarks. 
Compared to measured elliptic flow of electrons from D 
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meson decay in Au+Au collisions at y/s = 200 AGeV, 
the calculated values from the AMPT model, which is 
close to the D meson elliptic flow, are somewhat smaller 
even for an isotropic charm-light quark scattering cross 
section of 10 mb. Such a cross section is much larger 
than that given by the perturbative QCD estimate but 
is consistent with that based on novel resonant heavy- 
light quark interactions inside a strongly interact- 
ing quark-gluon plasma consisting of quasi-particles and 
colored bound states 0] • The study of charm flow in 
relativistic heavy ion collisions thus provides a promis- 
ing tool to explore the properties of produced partonic 
matter and its hadronization dynamics. 
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